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THEORETICAL AERODYNAMIC CHARACTERISTICS OF A FAMILY 
OF SLENDER WING -TAIL -BODY COMBINATIONS 
By Harvard Lomax and Paul F. Byrd 


SUMMARY 


The aerodynamic characteristics of an airplane configuration 
composed of a swept-hack wing and a triangular tail mounted on a 
cylindrical body are presented. For simplicity, the leading edge of 
the wing is considered to he straight and the trailing edge to he 
shaped so that the span-loading cirrve is flat between the fuselage and 
the wing-tip regions; the result is a nearly constant-chord avrept-hack 
wing. A method hy which other trailing-edge shapes can he studied is 
indicated. The analysis is based on the assun^ytion that the free-stream 
Mach number is near unity or that the configuration is slender. The 
calculations for the tail are made on the assumption that the vortex 
system trailing hack from the wing is either a sheet lying entirely in 
the plane of the flat tail surface or has completely "rolled up" into 
two point vortices that lie either in, above, or below the plane of the 
tail surface. 


INTRODUCTION 


The study of lifting surfaces flying at either subsonic or super- 
sonic speeds at small angles of attack has been reduced, by the well- 
known process of linearization, to the study of the equation 

(1-Mq^ cpxx Tyy + cpzz ~ 

where 9 is a perturbation velocity potential in a field having a 
uniform free-stream velocity Vq directed parallel to the x axis, 
and where Mo is the Mach niimber of the free stream. 

One basic simplification of equation (l) is brought about by 
neglecting velocity gradients along the span of the wing. If the wing 
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is lying in the z = 0. plane, this amounts to neglecting the tenn (pyy 
in equation (l), and results in the -well-kaown partial differential 
equation hy means of -which two-dimensional or section characteristics 
are studied. 

Another basic slniplification of equation (l) can be attained by 
neglecting the term (I-Mq^) Such a procedure is possible when the 

Mach number is close to 1 or the -wing plan form is so slender that 
velocity gradients in the free-stream direction are negligible in com- 
parison with the gradients in the y and z directions . Equation (l) 
has already been analyzed in these two connections in references 1 and 2 
for certain plan forms. The purpose of this report is to extend this 
theory, which has been named slender -wing theory, to include an entire 
airplane configuration. 

Results are presented for a nearly constant-chord, swept-back wing 
onoimted on a cylindrical, body having a triangular horizontal tail located 
eift of the -wing trailing edge. Both -wing and tail are flat surfaces, 
and the results are only those due to changes in the airplane angle of 
attack. 

A list of inportant symbols is given in appendix A. 


I - SWEPT-BACK WHIG ON A BODY OF REVOLUTION 

Partial Differential Equation, Boundary 
Conditions, and Form of the Solution 


Under the assumption that the free-stream Mach number is 1 or that 
the perturbation velocity gradient in the x direction is small, the 
partial differential, equation which must be satisfied for the solution 
of lifting surface problems can be -written 

^*yy + 9zz = 0 ( 2 ) 

Equation ( 2 ) is sinply Laplace’s equation in two dimensions, the 
variables representing lateral and vertical coordinates in a plane 
transverse to the direction of motion. 

The bo unda-r y conditions associated -with equation (2) eire given 
along a line in this transverse plane and specify that the fluid veloc- 
ity is everywhere tangential to the sinrface of the body. The problem is, 
of course, to find at other points in the plane the potential that 
satisfies equation ( 2 ) and fits these boimdary conditions. Of partlculeir 
interest is the streamwise coaponent of velocity along the surface of 
^ the wing and body since this is directly related to the loading thereon. 
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Solutions to equation (2) axe readily available. Two different 
analytic fonns of these solutions will be used in the fo3JLo-wlng analysis. 
One form is concerned with the use of the complex variable, the other 
with the use of Green's theorem and the inversion of a real, singular, 
integral equation. In general the procedure will be to use concepts 
associated with the complex variable to map the boundary conditions onto 
a slit along the real axis, then to solve the resulting problem by 
inverting an, integral equation, and finally, to use the complex variable 
again to extend such a solution out into space by the principle of ana- 
lytic continuation. 

Discussion of notation and transf onoations . - The first part of this 
report ttIII be devoted to the analysis of the configuration shown in 
sketch (a). The , f oUovring is a description of 
this configuration. Everywhere behind the 
leading -edge -fuselage Juncture the fuselage 
is a circular cylinder having a radius ro- 
Ahead of this Juncture the fuselage comes to 
a point, the manner being arbitral^. The 
wing is a flat plate without twist or camber 
movmted at zero incidence on the fuselage 
and the whole configuration is placed at a 
small angle of attack a with respect to 
the free-stream direction. The origin of 
the coordinate system is located at the 
wing apex. The leading edge is a straight 
line with slope, dy/dx, equal to m. It 
will be convenient at some places in the 
report, however, to use the expression 
y s(x) for the equation of the leading 
edge, hence, s(x) and mx are used inter- 
changeably. The trailing edge is repre- 
sented by the line y = t(x) and is, in 

general, not straight.^ The Tua-xiTmim semispan of the wing is denoted 
by So- The symbol to> as can be seen in the sketch, refers to the 
lateral distance from the x axis to the point at which the trailing 
edge intersects a line that is parallel to the y axis and passes 
throTigh the last outboard point of the leading edge. Finally, Cq is 
the chordwise distance from thp origin to the trailing-edge-fuselage 
Juncture . 



ilt was considered advisable at this time to consider only the rather 
particular configuration outlined. As the analysis progresses it 
i^ill be pointed out where the solution can be generalized to Include, 
for example, wings with twist emd camiber. 
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A second coordinate system will also be used in the succeeding 
development. Let tbe y,z plane be represented by the coD5>les: 
variable | , 

6 = y + iz == pe^® 
then Introduce the plane, 

li = + izi = Pie^^i 

so that the | plane maps onto the plane by means of the Joukowski 

transformation 

5 ^ ( 3 ) 

By means of such a transformation, the circle of radius Tq which 
represents a section of the fuselage in the I plane maps onto a portion 
of the real sixis in the plane (see sketch (b)); 

the part of the real axis which 
lies outside the circle^ in the £ 
plane maps into the remaining part 
of the real axis in the plane. 

As a consequence of eqxiation ( 3 )# 
the following relations hold for z^ 
equal to zero: 




w 



^The Joukowski transformation is double valued in that the regions 
inside and outside the circle p = r© both map onto the entire ^ 1 
plane. In this report only the field outside the circle is of interest. 
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and 


Further, 


2 



71 = 2to cos 0, 




®1 = ^ ti = t + -r^, ri - 2 to 


( 5 ) 


(6) 


From the basic theory underlying the use of complex variables in fliiid- 
flow theory, induced velocities in the two planes are related by the 
expression 


■iw = (vi-iwi) ^ 

d| 


from which, since in polar coordinates 


d| 


1 - 


cos 20 


+ i sin 20 


( 7 ) 


( 8 ) 


it follows that 


V = Vl 




w = w. 


V-r = 


1 — COS 20 


I 1 cos ii0 -V, ( 1 

J \PJ 

Vl cos 0+Wi sin 0 


sin 20 


sin 20 


1-1 ^ 


2 1 


+2wi 




) ( 9 ) 


sin 0 


Lastly, Laplace's equation must eilso be satisfied in the li plane, 
hence 




yiyi 


( 10 ) 
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Boundary conditions .- In this part of the report (part l), the 
effect of a cylindrical body mounted, on a nearly constEint- chord swept- 
back wing will be studied.. (Eeference 3 contains an analysis of the 
effect of a cylindrical body mounted, on a triangular wing, and reference 
presents results for a swept-back -vrlng with no body; both references 
use the assumptions of slender -wing theory.) The boundary conditions 
will be presented in the y,z space first and. can be -v/rltten 


(1) 

o 

II 

p = ro, O<0<2rt 



(ii) 

w = 0, 

z = 0, t^ < y“^< 8^ 


(u) 

(iii) 

<1 

II 

o 

w = Voa, p =“ , O<0<2 jt 

1 



Equations (ll) represent the conditions for a cylinder located at 
p = ro and two wing panels locate.! between ±t and ±s on the real 
axis, both cylinder and -brings being at rest in a free -stream which is 
moving with velocity w = YqCu at infinity. 


It follows from eq.uatlons (9) that these boundary conditions 
become, in the pleine, 

2 2 

(i) wi = 0, 0<yi < ri 

(ii) wi = 0, ti^ < yj^^ < si^ ' 

(iii) v^ = 0, Wi = Vott, =00, O<0<2« 


( 12 ) 


Equations (l2) represent the boundary conditions for three Wing panels 
along the real axis, aT 1 at rest with respect to the free stream moving 
with velocity w^^ = V^a at infinity in the transformed plane. 


It is more convenient to work with boundary conditions which vani sh 
at infinity, .Jiowever, so the final form of the conditions which must be 
satisfied is derived from equations (12) by subtracting the free-stream 
velocity Vga. There resiilts 


(i) 

wi = - Vott, 

0.<yi^<ri^ 

(ii) 

Wi = - Vott, 

t^^<y^^<s/ 

(iii) 

II 

II 

O 

pj^ =w , 0< 0< 2n 


(13) 


General solution .- The general solution to equation (lO) which 
gives the -rertical induced velocity w^ at a point in the li plane 
due to the Jump in the -ralue of the induced velocity Vj^ across the y^ 
axis can be written (see, e.g., reference 5) 
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(yi-yg)Avi(ya) 

(yi-ysf +Zi^ 


< 3^2 


m 


where is the variable of integration. Set equal to zero and 

there results the value of the vertical induced velocity on the y^ 
axis . Thus 


Wi(yi) = Wj^Cyi^o) 


AVj^(y3) 

2«j yi-ya 

-Si 


^2 


(15) 


Equation (15) is the form of the solution which will be used to 
analyze the problem previously outlined. It is apparent by reference to 
the boundary conditions listed as equations (13) that in equation (15) 
the value of Wi is the known quantity and Avi is the unknown. Hence, 
equation (15) is an Integral equation which must be inverted in order 
that the solution can be written. Such an inversion is not difficult if 
the value of Wi is known everywhere in the interval -si<yi<si. In 
the present case, however, there is a sub interval ri^<yi 2 <ti^ in 
which Wi is not specified, and further, in which Avi is not neces- 
sarily zero (due to the presence of a trailing vortex sheet). It will 
be shown in the subsequent development that the assumption that Avi 
is zero in this interval (i.e., no vortices are shed by the wing ahead 
of the interval) yields a nearly constant-chord, swept-back wing; with 
such a restriction the inversion can again be performed. 


Given the Inversion of equation (l 5 )^ it is possible to write both • 
Wi and V3_ for certain portions of the real eixis. All along this axis 
the functions w^ and Vj_ are, of cotirse, real. Hence, if 

f(^l) = Vi(yi,Zi)-lWi(y3^,Zi) (16a) 


then by analytic continuation 

f(gx) = v(ii,0)-iwi(gi,0) (16b) 

Therefore, the inversion of equation ( 15 )» together with equation (16), 
gives sufficient Infonnation to determine the induced velocities through- 
out space. 
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Particular solution for the nearly constant -chord vlng .- Adopt the 
notation equals Avi in the region of the axis representing 

the body or fuselage in the §3^ plane; Av^^o^ equads the vadue of Av^^ 
in the region of the sjcLb representing the space between the fuse- 

leige and the triing; and Av^^-w equals the value of Av^^ in the region of 
the yi axis representing the wing plan form. Then if AVj^o^ = 0 (the 
calculation of the traillng-edge shape corresponding to such a choice 
will be presented later), equation (15) becomes 


•w. 


i r 
J yi-7z 

-Si 


•ayg- 



yx-7s 





^•yiw 

yx-ya 


^2 


(17) 


Since the airplane is laterally symmetrical, the span loading is sym- 
metrical and A 9 i(yi) = A 9 i(-yi). Therefore, Avi has the property 
Avi(yi) = - Avi(-yi) . By means of this relation for Avi and the addi- 
tional change in notation 


Jli = 7x^ 
^2 “ ^2 


(18) 


equation (17) can be written 

_ 2 


w 


L(ni) - - tf 


%“^2 


Si 

-i/. 


%“ll 2 


(19) 


Equation (19) will now be inverted under the condition that 
Wi = - VoO' for 0<Tii<ri^ 2 

tional condition that (Av^^-^^) 


and for ti^ < t|i S i^, and under the addi- 


= (Avia) 


2=0, which amounts to 


Tl2=^x ^2 

assuming the Kutta condition along the wing trailing edge (see, e.g., 
reference 4 ) . This inversion is accomplished by a double application 
of the following solution (see appendix B); If 


f(%) - - i / 


AVi(T]3)dTl, 


( 20 a) 


then, iinder the condition that Av^^(a) = 0 

b 


= 2 /p! r /S a,, 

fi. ’ 


( 20 b) 
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Now write equation (19) in the form 


-Vott + 


tr 


%"'H2 


-tf 


%”^2 


and then, since Avj^-jj(o) = 0 by reasons of symmetry, apply 
equation (20h). For 0<Tij^<r^^, there results the expression 


Avib = - 2Voa 


Til 







dii2 

( 21 ) 


Substitute equation (2l) hack into equation (l9)^ reverse the order of 
integration and, for t3^^<Tii<si^, there results 


— j ^3- g(ng)<^^2 

2 k 111-^12 

ti^ 


(22) 


where 


g(il2) 



(23) 


Again apply equation (20h), this time to equation (22). 
g(Tli) can he shown to satisfy the relation 


In this way 


g(T]i) = - 2Voa 


Si^-Tll 


(24) 


and equating this expression to equation (23) gives 




ni(iii-'ti ) 

(si^-%)(Tli-ri^) 


f ti^<T)i<Si^ 


(25a) 
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A repetition of the above process yields’ 


Avib = - 27oa /- a 0<ni<ri^ 

(si -Tii)(ri -Til) 


(25b) 


The results given hy equations (25) can he extended to other points 

in space hy analytic Continuation. Equations (l6) indicate the necessary 

procedTore. Hence, since (vi) _ = Avi/2, 

z 


■Vi(yi,Z3^)-iW3^(yi,Zj^) = Voa 


— 

i + 



ax^(Si"-ti") 

i^-li=)(li^-ri^) 


(26) 


®Men the method is applied to a value of which has some given 

variation with tij_ there results 

' /t si^-Tl2)(ri2-Ti2) 


^ / T 2 / %(ti^-ili) 


o 


,1 Wi(t12) 

dTia . . - ^ 

■Tla V n2('ti -‘Ha) 


J Til- 


I 

ti" 


Si 


■WlCila) /(ila-ri^Ksi^-qe) 




dq 


Tl2(H2-‘tl^) 


, 0<Tii<ri" 


and 


^ f . 2 / Tn(ni--ti^) 


(si^-tIi)(t1i-^i ) 


Wi(ti 2) /(^i^~^g)(si^~^a) + 

_ J V ^ 2(^1 "^ 2 ) 


/ 




dq 


TI 1 -H 2 V ^2(^2"'^l ) 


, ti2<qi<si‘= 
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Equation (26) has several branch points so it is not uniquely defined 


without specifying the cut from -® 
the plane (see sketch (c)). 

In the upper half of the plane 

varies between 0 and it gind 
in the lower half, between 0 
and -It. Notice that when is 

at a very large distance from the 
origin in any direction the mag- 
nitude of the term on the right- 
hand side of equation ( 26 ) tends 
to zero, so that the boundary 
conditions at Infinity ajre satis- 
fied. It is evident that the 
other boundary conditions in 
equations (l3) are also satisfied. 


to si along the real ajcls in 





The Trailing Edge 


Special trailing-edge shape .- Equation (26) is a solution to 
Laplace’s equation and represents the flow around a wing and body. 
However, the plan form of the wing has not as yet been evaluated, 
although it has been fixed as that which makes the value of Avj. vanish 
in the region between the wing and the body. Since Avi is the 
gradient of Acpi in the yi direction, and further, since (Aq>i)m -n. 

(the value of A<T>i at the wing trailing edge) equals the total circu- 
lation about a given chordwlse section, this amounts to the same 

thing as assuming that there are no trailing vorticies between the wing 
and the body. It is a further consequence of such an assunrotion that 
the span loading ahead of this region is a constant for ro^ yo^ to^ . 
The configuration which will produce such a flow must now be determined. 
In particular, if the leading edge is taken to be a straight line, the 
equation for the trailing edge is unique and needs to be expressed. 

One of the simplest ways of finding the shape of the trailing edge 
is to find AT’t.e. f^rom equation ( 26 ) and solve for t as a function 
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of s for a fixed value of Tq and The constant represent- 
ing is the value of the jump in potential at the point P in 

sketch td). Here Acp is known (see, e.g., reference 3) since there is 

no gap to make its solution indeter- 
minate. 



Consider an- arbitrary section, as 
AA in sketch ( d) , downstream of the 
point P. The value of Av^^^ at such 
a section is given by equation (25a) 
and the solution for Api-w follows by 
definition and is 



(27) 


Equation (27) is an elliptic integral 
which can be easily reduced by means of 
the substitutions 


a 

ki =—2 2 

si‘=^-ri‘= 


> 


1 2 


4. 2 „ 2 
^^1 -^1 

si^-ri^ 


(28) 


and by using the Jacobian elliptic functions • defined, in this case, by 
(si^-ti^)sn^u = si^-yi^, cn^ = l-sn^, dn^ = 1-k^^sn^ 


to the form 


^*Piw “ 2700, 


2 2 

81 -ri 


E(ki,ti)-ki'^F(ki,ilri) 


(29) 


where 






(30) 


eind where the incomplete elliptic integrals E and F are defined in 
■the list of symbols (appendix A). Equation (29) reduces to the results 
given in references 3 ssid 4 when there is no gap or no body respectively. 
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At the trailing edge yi = t^, and equation (29) becomes 




2Voa ysi^-ri^ 



(31) 


where the elliptic integrals are now complete. Transform this to the 
physical plane, using equation (6), and set 


ko - 


y7s^t^-ro^(6^-t^ 


t(8^-ro^) 


, ko’ = yi-ko" 


then there results 


A9 


T.E. 


2VoO, 




(32) 


(33a) 


At the juncture of the fuselage and the wing trailing edge (the point P 
in sketch (d)), s equals mco and t equals Tq} so that 
equation (33a) reduces to 


(Acp)t.E. = 2Voomco 



(331) 


As was pointed out, the solution for the equation of the trailing edge 
can he obtained hy equating these two values of (A9)f]i . Hence 


or 


Set 




Eo-ko' 



t^ 


-mcot 


»• 


* 

< _ 

) ■ 

L J 

ko' 


Eo-ko I 


ro^ = 0 


V 


1 - 



Eo-ko 



(3h) 


G = mco 



Ik 
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and since s equals mx 



From equation ( 3 ^) tlie solution that gives the correct trailing-edge 
shape can he written 


t = 


G + 


ya^+i^ro^ 

~2 


(35) 


and from the definition of ko’ 


mx = 


.2 2 

t -ro + 




2ko't 


(36) 


If ko* and ro/mco are fixed, t/mco is determined from equation (35)j 
and a fixed ko‘ , ro/mco and t/mco determines x/cq from 
equation ( 36 ), Hence, it is relatively easy to ceilculate numerically 
the shape of the trailing edge. 

Sketch (e) shows the shape and position of the trailing edge when 
the wing leading e'dge is swept hack 45° and the radius of the fuselage 
is 31*6 percent of the extended root chord, cn* (A dimensionless 



coordinate system is chosen, 
however, so that the results 
can he used for various veilues 
of m, and Cq.) Shown also, 
for comparative purposes, is 
the position of the trailing 
edge when there is no fuselage - 
the condition in hoth cases 
helng, of course, that no 
vortices trail hack in the 
region directly behind the wing. 
Table 1 presents coordinates of 
the trailing edge' for several 
values of ro/mco* 


fej 
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Sketch (f) indicates the 
variation of the local chord 
along the span for Tq/to-Cq 
equal to 0.6, O.316, and 0. 

It is apparent that the 
effect of the body is to make 
the trailing edge more nearly 
that of a constant -chord wing. 
The asymptotic value of the 
wing chord is given by the 
equa.tion 



where 


c = 

“ x-» 00 \ ny 



rfj 


More conqplete results are 
given in table 2. 


Other trailing-edge shapes .- The procedure just presented can be 
generalized and used to calculate trailing-edge shapes corresponding to 
arbitrary span-loeiding curves. Sxjppose that the span loading in the 
transformed (or plane is represented by a power series in y^^ in 

the interval — Then the circulation in this interval can 

be OTitten 


r(y.) = ^ 

n=o 




and, hence, AVj_g^, the value of Av^ in the same interval, becomes 


Avia 


m 


n=o 


obnyi^"^ 


Equation (19) now takes the form 
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+ 


3 ^ 

2it 





T1 -Tl 
'l *2 





3 ^ 


The left-hand side of the latter equation varies with qx ^ given 
manner depending on the hn's in the expression for F. Hence, the 
equation can he considered as identical to equation (l9), the left side 
being regarded as an effective Wx in equation (19) . The analysis 
succeeding equation (19) can now be repeated in terms of the equivalent 
Wx. There results an equation for the trailing edge which depends on 
the b]x' s ^ 


By the process outlined, both the trailing-edge shape and the span- 
loading curve have been expressed in terms of m + 1 constants. By 
varying the number and magnitude of these constants, a large class of 
trailing-edge shapes can be obtained. 


The Wing Area 


Having foimd the shape of the trailing edge by the methods outlined 
in the preceding section, it is now possible to determine the area of 
the wing. Denote this area, region 1 in sketch (g), as Sx and the area 
of region 2, shOTm also in the sketch, as Sg. It is evident that the 
sum of these two areas is sin 5 )ly 
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If a dimeriBlonless system "based on the length Cq is adopted, one can 
write for the total area (i.e., hoth panels, see the shaded area in 
sketch (hi)) of the wing the equation 



where t/mco is given numerically as a function of x/cq in table 1 
(y/mco in the table representing t/mco) • Numerical relations between 
the parameters S/mco^, ro/mco# So/mco# and to/so are presented in 
figure 1. 


The area of a wing with another 
kind of tip shape can be readily 
evaluated once the particular tip 
shape is specified. For exanple, the 
area of the wing shown in sketch (h2) 
can be calcTilated by subtracting a 
rectangular area (given by the sum of 
the two triangular regions labeled 3 
in sketch (h2)) from the area of a 
sketch (hi) . 


Down^ra,sh Behind Wing 


The equation for the downwash 
behind the wing and in the z = 0 
plane follows Immediately from 
equation ( 26) . In the transformed 
plane the value of w^ is 

Wi = - Vou, 0<yi^<ri^;ti^<yi^si^ 

( 38 a) 


and 




(h) 


Wi = - Voa 


■ 

1 



ra^<yx^<ti^ 


(38b) 



18 


NACA ra 255^ 


In order to transform this value to the physical plane, care must he 
taken to go backwards through the boundary conditions in the proper 
order. Equations ( 38 ) represent the solution for the boundary conditions 
presented in equation (l3). To find the solution for the conditions 
given by equation (l2) & free stream Voa must be added. Thus, in 
mathematical notation. 


(■Wi) 


12 


(w, ) + VqU 

■^13 


where the subscripts 12 and 13 refer to the boundary conditions satis- 
fied. Finally, to find the downwash in the physical plane, the treinsfor- 
mations given as equations ( 9 ) aJ^d ( 5 ) must be employed and the free 
stream subtracted so that 


wCy) 


(wi)_ + V, 


13 


.a. 




V„a 


where (wt) now becomes 
•'■13 


(Wi)i3 " “ 0<y2<r^;t^<y2< 


and 



, r^< y^< t^ 


v(y) = - Voa,0<y^<r^jt^<y^< s^ 


(39a) 


and 


w(y) = - Voa 


1 - 


/s\ /( y^t^-ro^)(t^-y^T 

y^ J (s^y^-rp‘^)(s2-y2) J 


, r^<y2<t2 


(39b) 
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The accompanying sketch shows the variation of -w/Voa in the 
intervals for which it has been given. If no wing is a1;;tached to the 



body (or if the gap is very large) the fluid at the side of the body is 
moving upi^rard at a speed equal to that at ■which the body is mo'ving down- 
ward. The presence of the -wing restricts this motion and as the -wing 
panel approaches the body the air in the gap is forced more and more to 
move downward with the wing and body. The dotted lines in the sketch 
show the variation of -w/Voa if no body is present, that is, if ro 
equals zero. 


Chord-wise Load Distribution 


Loading on the wing .- The loading on the -wing can be calculated 
by means of the linearized equation for the loading coefficient. This 
equation can be written 


2 ^ _ 2 
Vn Vn 


w 


q. 


w 


(40a) 
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Xalin?ria.les s,, equation (^a) can .e 

modified slightly to read 

f ^ K 1 ^ ^ (UOb) 

V a V.^ 


w 


Bsi ds dx 


The value of 
thus 


can he obtained by differentiating equation (29); 


^lACp 


ITf 


Sb-l 


= 2Voa 




j^E(ki,^l)-kii2F(ki,Tlfi)^ + 


J 


^E(ki,ti) _ ,2 M^xM_ _ 2ki’F(ki,^l) ^ 

Ssi ^ 


which becomes 



The terms dk,/ds, and B^i/Ssi both involve dt^^/ds^^ which is 
proportional to the slope of the trailing edge in the transformed 5i 
plane. This latter derivative can be readily obtaii^d from equation (41; 
itself since the value of (Ap/q)^ and hence ^cp^v/dsx must be zero 
on the trailing edge, that is, where Jx equals t^ and ijr^ equals 1. 
This yields the relation 


'^'^1 ^ ^ 1^1 ( 1 ^ 2 ) 

dsx 
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by means of wMcb the Identities 


and 


_ Si / ^ 1-ti^ \ ' 

3sx si^-rx^ \ kx^i 





Sx^-rx^ \ Kx 




m 


can be written. Place equations (4l) and (k3) into equation (Uob), and 
there results for the loading 






i y(l-y/)(l-k/>|r/) _ 

i 

Transforming this to the I plane, one finds, finally, since 
dsx/ds = (s^-ro^)/s^ 


1 

m 


qa^ s^ 


E(ko,to) - - F(ko,+o) + 
, Ko 


s^(y^-rp^) /(y^-t^)(y^^-r^^) 

yt( s^-ro^) y ( s^-yS) ( B^y^-Vo^) 

where ko is defined by equation (33) and 




t /(B^-y^)(s^y^-r7) 
y y (B^-t^)(s^t^-ro^) 


(^5) 


In the special cases when there is .no body or when the wing is 
triangular, equation (Mf) agrees with the results presented in 
references (4) and (3), respectively. A discussion of the chordwise 
load distribution over the wing will be given at the end of this section. 
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Loading on the “body .- The variation of the load distribution over 
the body can be calculated in much the same vay as that over the wing. 
It is first necessary, therefore, to find the jun^) in potential between 
points directly opposed above and below the z = 0 plane. In the 
plane this difference follows Immediately from equation (25b) just as 
equation (27) 't-ra.s written for A<Piw. Hence, 


Acpib 



The first of the integrals in equation (46) has already been evaluated 1 

and the second can be reduced by means of the transformation 

(ti^-yg^) gn2^ = ri^-y^^ 

After some manipulation, equation (46) becomes 

Acpib = 2Voa rEi-E(ki,1f3)+ki*^F(ki,t3)-ki’®Ki 



where k^^ is defined by equation (28) and 



Using the equations (42) and (4 t), one can write for the loading on the 
body (after differentiation and simplification of equation (47)) 





NACA TN 255^ 


23 




2ko^(s^-ro2)t 


ro^-y^ 


[ ( 1^+ro^)2-l^y^t2] [ ( s2^ro^)2_l,y2323 

wtiere ko is defined by equation (33) and 


(49) 


ijr^ = 2t 


Again equation (49) agrees -vTltb 
previously knomi results in the 
limiting cases when the body 
vanishes or when the wing 
becomes triangular. 

Discussion of the chord- 
wise loading .- Equations ( 44 ) 
and (49) form the basic results 
of part I of this report. 

Graphs of the loading coeffi- 
cient for a wing alone and for 
a ^ring-body combination (ro/mco 
equals O.316) are shown in 
sketch ( jl) . The results for 
the case of zero body radius 
could have been obtained 
directly from reference 4, 

They are shown here for the 
purpose of a qualitative com- 
parison. Unfortunately, the 
load distributions on the two 
^•rings cannot be compared quan- 
tatively on the basis of equiv- 
alent plan forms since the 
trailing-edge shapes differ 
significantly. The variation 
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of the loadiBg along the center line of the hody is shown in sketch ( j2) . 



On the basis of the load distributions presented in references 3 
and 14-, the qualitative variation of loading shovm in sketch (jl) is 
obvious. That is, the loading falls steadily from its infinite peak at 
the leading edge to zero at the trailing edge. On sections which are 
cut by the Mach TOve from the trailing-edge fuselage juncture, the slope 
of the curve is discontinuous. 

The change in the load, distribution brought about by the presence 
of a wing tip is the same for a wing-body combination as for a wing 
alone. The behavior of the loading in the vicinity of a tip has a 
straightfoihrard explanation in terms of the trailing vortex sheet. 

Thus, if the wing is cut off along a line perpendicTilar to the free- 
stream directi cn, the vortices which were bound in the wing all turn and 
trail backwards with the same distribution in strength'^ as they had when 


^is assumes, of course, that the vortices have not begun to roll up to 
any significant extent. 
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crossing the last spanwise section of the wing 
( see sketch ( k) ) . Since the vertical induced 
velocity along the last spanwise section was made 
constant (by finding the appropriate solution to 
the integral equation ) , it must also be constant - 
and, in fact, the same constant - everywhere in 
the vortex wake. Hence, if a flat stirface having 
the same angle of attack as the wing is inserted 
anywhere in the wake it will in no way disturb the 
flow and consequently there will be no loading on 
such a surface ( just as there is no loading on the 
vortex wake itself). The loading is zero, therefore, 
on the tip regions marked 2 in sketch ( k) ; the load- 
ing in the regions marked 1 being given, of course, 
by equation ( 44) . 

It is interesting to see how the distribution 
and magnitude of the loading given by this 
( slender wing) theory con 5 )are with linearized 
theory results at some Mach number other than 1. 

The differences caused by considering Mach numbers 
other than 1 depend, of coiurse, on whether or not 
the new Mach number is subsonic or supersonic. 

This discussion must be limited to a coniparlson 
with supersonic Mach numbers only, since theo- 
retical chordwise load distributions over swept-!- 
back wings flying at high subsonic speeds are not 
available. The change in the loading brought about by 
increasing the speed can be divided into two parts: 
one, a change caused by the rotation of the Mach 
lines which form the boundaries of the various 
regions in each of -tThich the shape of the loading 
cu 2 Tve rakes widely different forms; and the other, 
a change in the magnitude of the loading within 
each of these regions. 



(k) 


Sketch (Z) indicates these effects. Thus, on the wing flying at 
supersonic speeds the sharp drop in loading occurring at a critical Mach 
line moves farther back along the chord from point b to point a in 
sections AA and BB shoTO in the sketch. This causes a considerably 
higher value of the loading for the supersonic i-rlng® in regions 1 and 2. 

A similar effect occurs on the body traveling at a supersonic speed where 
now, however, the traces of the Mach lines are no longer straight but. 


^Solutions showing the effect of crossing critical Mach lines on a 
swept-back supersonic wing are given in references 6 and 7» 
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I. — Supersonic 



0 c 


(!) 


due to the curvature of the body, form 
helices. Region 3 in "the supersonic case 
would be a region of zero loading and region 4 
would be a region of high loading relative 
to the sonic value. 


The relative magnitude of Ap/q.a vriLthin 
the various areas bounded by the t'ring edges 
and the pertinent Mach lines changes as the 
reference section moves outboard along the 
wing span. Along inboard sections ahead of 
region 1 in the sketch (i.e., ahead of the 
sonic Mach line from" the trailing edge root) 
the loading on the sonic wing is higher than 
that on the supersonic wing. It is well 
known, for example, that in the case of a 
triangular wing mthout body, slender -wing 
theoi*y gives a loading E times rhe loading 
obtained at a supersonic Mach number (where 
E is the complete elliptic Integral of the 


second kind with modulus i/l-m^P^ and is given 

closely by • 1 + i m^P^ (in — - ^ for 

2 \ mp 2 / 

small values of mp) . On the other hand, 
along sections farther outboard, the magni- 
tude of Ap/qa on the sonic mng must become 
lower than that on the supersonic wing. This 
follows immediately from simple sweep theory, 
since the component of velocity noirmal to the 
leading edge is closer to the speed of sound 
than xhat for the sonic wing. In fact, it is 
easy to show that at distances far enough 
outboard so that simple sweep theory applies® 
the supersonic wing has a loading 


i - 


times that obtained from 


slender wing theory. 


By an application of the above consider- 
ations, it is possible to obtain an estimate 
of the absolute value of the loading on a 
wing-body combination ax supersonic Mach 


® Sketch (jl) indicates the manner in -which the loading approaches that 
given by simp)le sweep theory as "the reference s-tation moves outboard. 
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nTjmbers. Another manner in vhich. the results of slender wing theory can 
he extended to Mach numbers other than 1 (or to pl^ forms which are not 
sufficiently slender) is to form the ratio of the resulting values for 
the Tfing plus body to those for the wing alone and apply this ratio to 
solutions for the same wing or body at the required flight Mach numbers 
(or slenderness factor). As has already been mentioned, the formation, 
of such a ratio for the load distribution is not possible from the solu- 
tions presented herein since the wing trailing edges change to a certain 
extent with the addition of the body. It is reasonable to expect, how- 
ever, that a ratio of the Integrated loading characteristics (i.e., lift, 
drag, and pitching moment) formed by dividing the result for a wing -body 
combination by those for a wing alone -sfill be useful in estimating the 
Interference effects even if the xvLng trailing edges differ slightly. 


Aerodynamic Characteristics 


The results developed in the preceding section can now be converted 
into forms which represent the aerodynamic 'characteristics of the wing 
and body. Hence, the following will present the span loading, average 
chord loading, lift, drag, pitching moment, and center of pressure for 
the wing-body combination. 

Span loading. - The development of 'the span loading on the wing and 
body will be considered separa-tely. First the span loading on the wing 
can easily be determined from the value of AT given in a preceding 
section. Thus 

// . 1 ^^- 

wing plan form 

and since AT at the leading edge is zero 

span 

where (AT)ip^j;^ is the -value of AT on the trailing edge of the wing. 
Since (AT)ij g eilso represents the total circulation about the wing 
chord, there resTilts for the circulation developed by the wing and 

the total wing lift 
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The equation for can now be 

determined from equation (29). 

Between y = ro and y = to (see 
sketch (m)) the value of is a 

constant (this being the condition 
by which the shape of the trailing 
edge was determined) . Between y = to 
and y = So, ^w is given by the veilue 
of A^> along the section AA since 
there is no loading between this 
section and the trailing edge. Hence 

= 2Voacom f ^o!$y^'to 

(52a) 

and 


B ^~r ^ \ 2 

2Voa(^ -° -— J |^E(ko,io)-ko‘ F( 5 ^o,'I'o) J » to<y<so ( 52 b) 

where ko and '^''o are defined in the table of synibols. 

Some care must be taken in order' to find the span loading on the 
body. Since we are concerned here ^dLth the loading developed behind the 
wlng-leadlng-edge fuselage juncture, it is necessary to subtract the 
value of A 9 b ‘this station, shown as s'tation AA in sketch (n), from 
Aq>b at station BB also shown in sketch (n) . For the to-bal span loading, 
then, it ■will be necessary to add to ■fchis value the load accumulated 
on the nose of the body. Denote by (rt,)^ the increment of circulation 
developed by the nose of 'the body and by (r^ )^ the increment of circu- 
lation developed behind the wing-lead 1 ng-edge fuselage juncture, "that 
is, be-tween s'fcations AA and BB."^ Hence, 


"^Slender ■wing "theory gives zero loading behind station BB as long as 
the trailing vortex pattern does not vary. 
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The VBJ-ue of (Acp)l^j;.> sketch (n) , 
caxi he obtained from eqiiatlon (^7) hy 
setting ti and si equal to ri. 

By transformation of the result into 
the physical plane, one obtains 





The spEin loading on the body is then 
given by 




Eo-E(ko,T)f2)+ko’^F(ko,t2)-ko’^Ko 


+ 


2 


So 


/[(so^+ro^) -^8o^®](ro^-y®) 




2„2 


- 2 



and ijra 


where kQ 


are given in the table of symbols. 
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Voocom 



table 2. Typical resiilts are sbo^wn for a 
and 0.316 in sketch, (p). 


Sketch (o) shovrs the 
variation of the span loading 
over the wing and body for a 
body radius factor, ro/mco. 
equal to O .316 and a wing 
semlspan factor, so/mco 
equal, to 1.7* 

Section lift .- The wing- 
section lift coefficient can 
be calculated readily by 
dividing the section chord 
into the value of the span, 
loading at the same span sta- 
tion. Along sections not 
influenced by the tip cut-off 
this is especially simple 
since the span loading is 
constant. The VEilue of the 
section lift curve follows 
immediately, therefore, from 
body radius factor equal to 0 



(P) 
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Section drag.- The value of the section drag can be written 


cd. = acj + (c^)^ 


( 55 ) 


where (cd.)g represents the suction force at the section leading edge. 
The magnitude of (cd.)g can be evaluated (see, e.g., reference 8) by the 
equation 


(cd)g 


1 . ® 

qc dy 


(56) 


where dP/dy is the suction force in the free-stream direction per unit 
length normal to the free stream, and c is the local chord. 


Define a new set of coordinates, as sho™ in sketch (q), such that 
y^ lies along aud lies perpendlciilar to the leading edge of one 

wing panel. Then if 


r - It'll a 


un(xn^yn) 


Vr 




(57) 

the suction-force con^ionent F (positive 
in the positive x^ direction) in the 
free-stream direction is given by the 
eqviation 

S = -^G^(x,y) ( 58 ) 

<ay Pn 



Now by differentiating equation (29) with respect to x and dividing 
by 2 (to convert Au into u) , there resiolts 


u = mVo 




Ko yt(s2-ro^) y (s^-y^)(s^‘'-ro^) 


y/mco > 1 


k 
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and similarily v "becomes 


V = - "Voa 


y +ro 


©/ 


(y2t2_ro^)(t2-y2) 


y^ \t/ V (s%2_p^4)(gS_y2) 

Since the normal con^ionent un is given "by the equation 


, y/mco> 1 


Un = 


mu-v 




there results 


Ur, = m' 


Voa ^sS+rpg ^ 


■“ y 


l+m*^ 


E(ko 


r m^(s^+ro^)(y^-ro®) 


B^-ro^ 


By means of equation ( 57 )> G- 
Pn = l/r/ 1 +m^ , 


^ -* yl^l-Hn^ J ( s^-y^) ( s^-ro^) 

y J 

can now he calculated. Hence, since 


^ _ g(y^+ro^) / (y^-t^) (y^^-r~^ 

yt(l+m2)^/^ y . 2y(y^ro^) 

Finally, therefore, the suction force can he written 


dF _ a^(y^+ro^ I" (y^-t^Cy^^-ro'^) ‘ 

ay' L. 2y(y^-ro^) J 


(59) 


and, hy using equations (55) and (56), "the section drag coefficient 
he written 


can 
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Cd 


cm 


rt(y/mco) 

(c/co) 


3^0 


2-1 



y/mco>l 


(60a) 


In the region "v^here the leading-edge suction force cannot he affected hy 
the traillng-edge shape, equation (6o) reduces to the simpler form 


, r6/mco<y/mco<l (60b) 

The variation of the section drag coefficieht is shown in sketch (r) 
for two wings: one without a body, and the other with a body radius 
factor equal to O.316. 


Cd ^ £1 _ ^(yAco) 

cAn rrm ( L 






3^ 


MCA TIT 255^ 


Total lift . - The lift on the -VTing, the "body, and the conibination 
can no\r he evaluated hy means of equations f5l) and (53) since the 
expressions for f are given hy equations (52) and (5^). The integra- 
tion required is somewhat involved algehraic al l y , hut the final result 
can again he expressed in terms of elliptic integrals. Thus, defining 
AQ(k,T|r) hy 


Ao(k,T{r) = I KE*(k,iif) + EF'(k,T|f) - KF’ (k,T^ ) 


(61) 


( tabular values for Aq can he found in reference 9)^ the total lift 
carried hy the wing is 


qa 


2 / tp^+ro^ 

Sp \ tp 



-hXr 


‘ 

■ Ep-kp* 





It 



Ap(k2,^4)+2it 



(62)' 


Equation (62) agrees with the results presented in references 3 ^ 

when tp = rp (the case of a triangular wing on a body) and rp = 0 
(the case of no body), respectively. Men Sp eq uals tp, that is, 

■vdien there is no wing, Ly equals 0. 

The lift on the body will he computed in two parts just as was the 
span loading on the wing: the lift on the portion of the body behind the 
■wing -leading-edge fuselage Juncture (Lb)i, and the lift on the nose of 
the body (L'b)^. It is a well-known res'ult of Munk’s airship theory that 
the lift on the pointed nose is Just 


(L-b)p 

qa 


2itrp^ 


( 63 ) 


anfl is independent of the shape of the nose.® The ‘value of the lift in 
the -vicinity of the wing follows from the Integration of equation (5^) 
according to equation (53). The total lift on the body can then be 
written 


®In this report, it is assumed that the nose is always ahead of the wing, 
that is, the portion of the body on which the -wing is mounted is every- 
where a circular cylinder. 
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qa 


qa 


So 


Ifro { Eo-ko’ Kq 


to^+ro^ 


Eg-ko ^2 


k 



2itro^ 


(64) 


Setting to = Tq, one findB the. result given in reference 3 for a trian- 
gular wing mounted on the cylindrical portion of a pointed "body. 

If to = Tq = Bq, the -vring disappears and equation (64) reduces® to 
equation (63)« If r^ = 0, Lj^ reduces to zero. 

Finally, the Siam of equations (62) and (64) gives for the total 
lift of the wing-body combination, including the nose of the body. 


L. = + tJ) 

q.a. V so^ -t^^ °y 

(65) 

Sketches (s) and (t) show 
the total lift on various wing- 
body combinations together with 
its division into the component 
parts carried separately by the 
iTing and body. Various lift 
coefficients, depending on the 
choice of the reference area, 
can be fonmilated by means of 
the area-span relationship 
given in figure 1. 


Total drag .- In general, 
the voiiaex draig can be calcu- 
lated by finding the momentum 
transport through a plane perpen- 
dicular to the X axis and 
located infinitely fax behind 
the airplane. In slender wing 
theoiy the calculation of the 
total drag is simplified in two 



®Note that AQ(k,l) = 1. 
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ways: first, the vortex drag becomes the 
total drag (neglecting, of coxirse, vis- 
cosity), and. second, in the calculation of 
this drag the reference plane can he located 
immediately behind the airplane since the 
flow there is the same as it is infinitely 
far hack. 

Hence, a momentum balance gives for 
the drag 



Body 


D = - Po y' WoAfPdy+poTo 




(9) d0 

r=To 3 ^=To 


(66a) 


where w is the value of the vertical 
induced, velbcity behind the wing in the 
2=0 plane. It is more cqnvenient to 
perform this integration in the plane. 

Equation (66a) can be put in the form 


D = 


3 X* 

— r * AtPiSyi-Pt, r i Acpj^ayi 

^ \ V I 


(66b) 


where w and A<Pi are given in the following. 

For ri^<yi^<tx^, that is, between the body and the trailing edge, 
it is seen from equations ( 39 ti) and ( 31 ) that 


w = - Vcjo, 


1 - 2yi 




1 / ti^-yi^ 

-I J Bi^-yi^ 


Aq>i = 27oa 


(Ex“kx 



> ( 67 a) 
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For that is, on the wing, it is seen from equations ( 39 a) 

and (29) that 


w = 


Acpi = 27oO- - ki* F(ki,ti) j 


(67b) 


Finadly, for that is, on the body, equations (38a) and (^ 7 ) 

give 


w = -Ycp^ 

AtPi = 2 Voa ysi 2 _r ~2 


Ex“®(^i^^a) ^1* ~ 


V (67c) 


k 3 _’^K^ 


+ 2 VoO- 


(si^~yi^)(^i^~yi^) 

t ^-v 2 
-y^_ 


The substitution of equations (67a), (67b), and (67c) into 
equation ( 66 b) yields after Integration 


D ^ 

qa^ 2 qa 


- Ksi2-r^2)(E3^-k/ K3^)(E;^’-ki^2K^’) 


(68) 


where L/qa is given in equation (65)* In the | plane equation ( 68 ) 
can be written in the dimensionless form 


D 


kqa^SQ^ 8qaso‘ 


/ \ 2 

- ( ~l —J (Eo-ko’\)(Eo’-ko%’) 

k So^ 


(69) 


which for rg =■ 0 agrees with the res\ilts of reference k. Equation (69) 
adso checks with the result obtained for the drag by the method, pre- 
sented in the preceding section on section drag, based on the calculation 
of the suction force along the wing leading edge. 
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Sketch (u) shows the 
totaJ. drag on various I'ring- 
"body coiribi nations . 

Chord, loading .- In 
ord.er to find, the center of 
pressure and. pitching moment, 
it is convenient to find, 
first the chord, loading, 
which we will define as the 
value of /(Ap/qa)dy where 
the integration is carried 
over the wing and body. The 
chord loading can also be 
obtained by evaluating the 
expression d(L/q.cc)/dx 
since the latter term is 
equal to /(Ap/qa)dy. As 
a check, both methods were 
used to derive the following 
expressions- 

For the part of the 
chord loading contributed 
by the wing it can be shown 
that 




4(so^-ro^) 


2soTo 
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(70b) 


For the part contributed by the body behind the leading-edge fuseleige 
Juncture it can be shown that 


^ f S ^ Bin -2£or^- ^ /,,<x<c 

J \<1% ^ q-a so3 So^+roS " 


3 r f 
j \ ^Ay 


dy = 


d ^(sQ^-rp^) 


dx qa 


m 


(71a) 



\ AD(ka,t4) 




(71b) 


The total chord loa ding can be obtained by combining equations (70) 
and (71) • There results the e^^jressions 

for ro/m<x<Co 



and for Cq < x 



( 72 a) 


(72b) 
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Sketch (v) shows the variation of the chord loading with x 
for ro/mco equal to O.316. 



Center of pressure.- The results of the last section can he used to 
determine the center of pressure x<,^p^. The value of x^.p, is given 
hy the equation 


^c.p. 





(73) 


which excludes the loading on the nose. By means of equation (73) 
sketch (w) was constructed. 
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II - ADDITION OF A HORIZONTAL TAIL 


It is possible to use the calculations given in the first part of 
this report to find the forces and moments induced on a horizontal tail 
by the presence of the wing and body. The same assunptions that were 
used for the solution of the load distribution over the wing and body 
will be made here. Hence, the results idLU be principally valid for air- 
planes having highly swept -vrings and tails or flying at Mach numbers 
approaching 1. 

In addition to the basic assun 5 )tions by which slender, wing theory 
is defined, however, some additional assumptions must be made concerning 
the behavior of* the vortex sheet trailing behind the wing and passing by 
the tail. Actually these trailing vortices provide the only means by 
which the wings can signal their presence to the tall, and except for 
them the slender wing theory analysis of the tail effectiveness would be 
identical to that described in part I for the wing. Only two types of 
trailing vortex patterns will be investigated. One composed of a flat 
vortex sheet situated entirely in the z = 0 plane (the plane of the 
wing) , and the other co^i^osed of two coii 5 )letely rolled up point vortices 




MCA ra 255!)- 


k2 


situated symmetrically witli respect to the y = 0 plane, located a 
distance h above the z = 0 plane and a distance a from plane of 
symmetry. These patterns represent the two extremes of the actual phys- 
ical behavior of a trailing vortex wake. It is to be expected that the 
sheet is more representative of the true wake when the tail is located 
only a short disteince behind the wing. On the other h a n d, the two point 
vortices should be valid for tails located a large distance b ehin d the 
wing. An indication of the magnitude of the distances at which the 1 wo 
assungitions are accurate can be obtained in reference 10 . 


Method of Solution and Boundary Conditions 


The partial differential equation that governs the flow in the 
vicinity of the tail is, of course, identical to the one studied in the 
first part of this report, namely, Laplace’s equation applied to a yz 
plane (equation ( 2 )). In fact, the general discussion of boundary con- 
ditions and forms of solution given in part I still applies here. Hence, 
the Joukowski transfomnatlon can again be used, the | plane having the 
same relation to 'the plane as before and the integral relationship 

given as equation (15) still applying. 

The only mathematical difference between the study of the wing and 
tail can be seen at once in the application of equation (15) • In the 
case of a flat wing, the vertical induced velocity w in equation (15) 
was known to be a constant over the region occupied by the wing plan 
form. In the case of a flat horizontal tail, on the other hand, the 
value of w over the region occupied by the tail plan form is composed 
of two parts; one, the constant value fixed by the inclination- of the 
surface to the free stream, and the other, a distribution that is Just 
equal and opposite to the vertical velocity Induced over the region by 
the vortices trailing from the wing. Effectively, therefore, the analy- 
sis of a horizontal tail is the same as that for a wing with a given 
variation of t^rf-st and camber. 

The additional notation necessary for the description of the 
pertinent tail parameters is shoam in sketch (x) . The distance from 
the X axis to the tail leading edge is represented by a, and the 
slope of the tail leading edge is designated by d. In this report, 
only triangular tail shapes will be considered; however, more complicated 
shapes could be analyzed by the method presented. 
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Solution for Trailing Vortex Sheet 


Since the vortex sheet from the 'vring is assumed to lie entirely in 
the = 0 plane, and since the outer extremities of this sheet are 

at ±Sj^ (see sketch (y)), the study of this case can commence vith the 
inversion of the integr^ equation" (15) of part I. 


a i2< y t ^ 2 ^ 


Wi(yi) 


— p 

J yi-Va 

-Si 


dya 


( 74 ) 


the value of Avi(y2) is given hy equation (25a) and 


For 

for 


AVi(y2) = 0 


( 75 ) 
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Transformed plane 


Wing 



vortex wake 


7 


S. 


y, 


(y) 


Substitute these tiro values of Av^ into equation (7^) and apply the 
boundary condition that wi equals -VoO' da the interval 0 <yi 2 <Cj^ 2 . 
Then, assuming ti^>cfi^, that is, the vortex sheet from the wing does 
not cross the tail^® (the condition shown in sketch (y)); there is 
obtained, after inversion (see appendix B) and some manipulation, the 
val.ue of Avj^ on the tail . Thus 

for 0<yi^<cTi^ and for t^^^ > 0 ^^^ 


Avi(yi) 


j 2 2 

-yi 



/ (Ti-ti^)(T]-ai^) 

y (s3^^-Ti)(Ti-r^^) 


(76) 


This solution for Av^ can now be used, to determine the aerodynamic 
characteristics of the tail in the presence of the vortex sheet trailing 
back from the wing. 


^°This assun 5 )tion applies to a.1 1 subsequent a na l y sis of the tail and 
vortex sheet combinations. 
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Spem loading .- The spanwlse variation of circulation generated hy 
the tail surface is given by the expression 


A9^=^^Av(ys) ay^ (77) 

where Av is given in the plane hy equation ( 76 ). On the portion 

of the plane that is covered hy the tail surface (i,e., 

for ri^<yi 2 <(j^ 2 ) this yields 


Aqp^^ = 27oa 




arc tan 




dT] 


(78) 


To detenaine the span loading on the hody it is necessary to subtract 
the value of A<Pit a,t the tall-leading-edge fuselage juncture. This 
value is obtained from equation ( 78 ) by placing load- 

ing on the body (i.e., for 0 ^y^^^ <rj^^) is then given by the formula 



Equations ( 78 ) and (79) have been transformed to the | plane and the 
results are shown in figure 2(a). In this and in all following numeri- 
cal examples the wing will be fixed as the special type studied in 
part I having the meas\ires Bq/mOq = 1.7? "to/^'^o = 1»091 and 

ro/mco = 0 . 316 . 
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Chordwlse load distribution on the tall .- The distribution of load 
over the tail surface can be calculated from the equation 


/Aj)\ _ 2 _ 2 clAfPit 

YqO. dx 


(80) 


where the value of determined from equation (78), is given 

by the expression 


8AcPit 

Bcf, 


2YcfUSj^ p 

It y ffi2-yi2 J ^ y( ( tl-r^e) ( n-di^) ( Ti-t^^ 


2VoOTi 


2 - 




2) 


By means of the substitutions 


(81) 


,, 2_ _2„ _ (si^-(Ji^)(Ti-ti2) 

it ^ — 1 ^ Bli U — " ' " ■ ■ * 


(si^-^i^)(ti^-ri^) 

equation (8l) can be reduced to 




( 82 ) 


It V 


where 


= 




Hence, the -loeiding on the tail can be written in the closed form 






and since 


(83) 


(84) 


(85) 


££i ^ 
dx 
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the final expression obtained by transforming equation (85) to the 
5 plane is 






( 86 ) 


where k4 and '^s are the transformed values of kg and \|rg and are 
given in the table of symbols. 


Total lift on the tail .- The total lift on the tall can be eveiluated 
by use of the equation 


h 

qq, 



which becomes in the li plane. 


(87) 


Lt 

2qa 


^ / A^(y,)dy, + ^ f 


yiACfa.(yi)dyi 
V yi^-ri^ 


(88) 


The value of A<Pi is given by the equations (78) and ( 79 )* Substitute 
these ejqiressions into equation (88) and, after integration, there 
results 


^ = n(ti2-ri2-si2+0i2) Ao(kg,+g) + 


2 y(si2.ai2)(ti2-ri2) fEg 

L V St -St / 


(89) 


In the I plane this can be written 


2qa 


= — rt 


( So^-'to^) ( So^to2_ro4) ( 


2_, 


S 2t 2 
O O 




Ao(k 4 >^s) + 


2(to^-ro^) 

"to^O 


\/ ( So^CTo^) ( So%o^^o'*) ^4 " 1^4^ 


(90) 


A plot of this equation is shown in figure 3 (a) for ro/mco = O.316. 
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Dreig.- The drag of the tail can he calculated from the equation 


Dt = aJLt + F^. (91) 

where is given hy equation ( 90 ) F-t is the suction force at 

the leading edge. In a manner similar to that given in the first part 
of this report under the subheading "Section toag," this force is 
obtained from the equation 


Ft 



As before 


lim o un(xn>yn) 

^ Pn ^ 

Xn-^ o Vq 




(92) 


(93) 


the values of u^ and x^ being, respectively, the normal velocity to 
and the normal distance from the leading edge. Substituting these values 
into equation (93) gives 


G 


2 






(94) 


and the expression for the suction force can be written 


Ft = - VoVprt jT (~^) Ao^(k4,+s)'3y 
^o 

Hence, the total drag becomes 


(95) 


^ ^ Lt _ 

a^q. oq. 




^o^(k4> ^s)'ay 


This equation is plotted in figure 4(a) for r^/rncQ = 0.316. 


(96) 


Chord loading .- A closed formula for the loading can be obtained by 
carrying out the integration / (Ap/qa)tdy over the tail or by evaluating 
*3-(lit/<ia.)/dx. However, since iiie term Ao(k 4 ,+s) ^ expression 
for Ap/qa does not involve y, it is easier to evaluate /(^/qo.)t^» 
Thus, using equation ( 85 ) ^ the chord loading in the Si plane is given by 
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2 dy = rY'^Vi + 
j VWt VWt J -/7IW“ 


= W3_A^(k^,t9) 


dtJ^ 

dx 


(97a) 


In the I plane this becomes 

" ' w(2^‘) A,(k„t 3 ) (97b) 

o 

(This checks with the result found by differentiating Lt/qa with 
respect to x) . 


Center of pressure .- The center of pressure can be calculated by 
means of tha fn- rrmilp 


I 


cfo/Vi 


^.p. “ 


^o/Vi 


<3Lt , 

X — i dx 
dx 


Lt 


(98) 


where dLt/ dx is the resxilt just obtained in the previous section 
and Lt is given by equation (90). Placing these values into 
equation ( 98 ) yields 


X 


c.p. 



Ao(k4,’l^8) 


(99) 


A graph of this result is shown in figure 5(a) for ro/mcg = O. 316 . 


Solution for Rolled Up Vortices 


As was. pointed out in the preceding section where the method of 
solution was discussed, the manner in which the present problem will be 
attacked is as foUbws; First, the velocities induced at the surface of 
the wing and body by two point vortices located somewhere in space will 
be calcTilated; second, a solution will be formulated (by methods identi- 
cal to those uBed in part I of this report) that will just cancel the 
vortex induced velocity con^ionent normal to the surface of the wing or 
body; third, an additional solution will be formulated that will fit the 
boimdary conditions prescribed for the tail surface (in this r^ort only 
a flat-plate tall surface will be considered) . 
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As usTial, it is simpler to work witti the transformed plane 

than with the physical | plane. Hence, again the Joukowski tremsfor- 
mation will he applied to the field equation and boundary conditions at 
the outset of the problem. See sketch (z). 


i iz Physical plane 



The velocity potential 
at a point (yi,zi) in the 
1 3^ plane Induced by a pair of 
point vortices located at 
yi = ±ai and zi = hj. is 
given by the equation 

cn = ^2 (zi-hi)dyg 

^ J (yi-y2)^+(zi-iii)^ 

( 100 ) 


, Transformed 


plane 


Vortex 


(z) 



location 


7 


where F-w is the strength of 
the circulation carried by the 
wing panels and trailing back 
from the wing tips. The value 
of F^/VoCtmco that corresponds 
to the swept-back wings studied 
in part I of this report is 
given by equation (32a), thus 



The values of a^^ and hj_ depend, of course, on the span of the wing. 

In order to con5)are the following results for the rolled up vortices with 
those obtained for the sheet vortices in the preceding section, the 
numerical results presented in the succeeding examples will be for 
So/mco = 1 . 7 , to/mco = I.09I, and ro/mco = O.316. This particular 
choice of parameters fixes the span-loading curve for the wing to be that 
shown in sketch (o). A reasonable choice for the value of a can be 
calculated by replacing the figure in sketch (o) bounded by the lines 
y = Tq , r = 0 and the curve for F /VQOtmcQ by a rectangle with the same 
height and area. The value of a is then given by the sum of the base 
length of this -rectangle and the quantity ro/rncQ. This procedure was 
carried out and the result a/mcg = 1 « 5^5 'was obtained. In order to 
obtain a more conplete picture of the effect of the point vortices on 
the aerodynamic characteristics of the wing, four different locations 
were chosen for the positions of the vortices, two in the z = 0 plane 
for values of a/mcQ equal to 1*5^5 1.3, aad "two at a height 

h/mco = 0.3 above the z = 0 plane for the same two values of a/mco- 
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From equation (lOO) it- can easily "be shown that the value of the 
vertical induced velocity in the = 0 plane is 



yx-ax 

yi+ax 

■^^-i4x=o " 2it 

Lhx^+(yx-ai)2 

bi^+(yi+ax)2 - 


The particular inversion of equation (l5) that fits the present boundary 
conditions (see appendix B) can he written 


A'i'i(yi) = — ■ ^ ^( 3 ^ y dye (102) 

Place the value of Wv given by equation (lOl) into equation (l02) and 
add the condition that the tail be a flat-plate lifting surface at an 
angle of attack a, and there restilts 


AVi(yi) = 


-2Voayi 


Tw 


y ai2_y^2 2«y di^-yi' 




__ — + 


yi+^x 


yi+bi 


where 


y hi^« 




ybi^-cri^ 


yi-^] 


yi-hj 


(103) 


hi = ax + ihx 
hi = ai - ihi 


(104) 


and where the radicals are defined uniquely if the complex plane is cut 

along the real ajcis between -<» and Ox. (For example, can 

be set equal to where <Px must lie between -jt and it.) 

Although the above expression for Av can be put in real terms by 
applying the transformations 

^1 == cosh cos o)x 


|hx I = Ux slhh 7x sin Wx 


( 105 ) 


it is easier in deriving subsequent quantities from Avx to use 
equation (l03) first and to make the transformation afterwards. 



Span loading. - The span loading can he determined from the relation 


ro 


where Av is given hy equation (l03)» There final ly results 
/ _ 2 

Acpit =• 2Voa A/hi^-y^a arc tan 

rt 

*J [ ( slnh^73^-cos^W3^)+y]^®]^+4a3^'®(ax^-yi^) sinh^yj^cos^a^ -<Ji^(Blnh®7j_-cos^(jj]^)-yi^ 


(106) 


2n 


1 sinh COB 


yi^< 


( 107 a) 


and on the body 


A9it =■ 2Vott ^ arc tan 


y -cos^wj^ ) +y 1 ® +^i^( cTi^-yi®) slnh^/^cosSci - o,^( 0lnh£7i-coB2w3^)-y^2 2 Tw 

— — ' — + ~ arc tan 

' 2di 1 ^-yi® slnh 7i cos 0)1 


slnh^7£-co82o)2)+yi^J ^+^ri^(ri^-yi^) slnh^72Cos^2 -ri^(Blnh?7i;-cos^(»)^)-yi2 


2i^ <J slnh 7^ COB co. 


, Ti=<r3_2 


(107b) 
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The terms sinh 7^, slnh /g, sin o)i, and sin Ug are all defined in the 
table of symbols. Equations (107a; and (lOTb) can, of course, be trans- 
formed to the physical plane by transforming each symbol therein from 
the system to the | system. 

The variation of the tall span loading, as given by equations (l07a) 
and (107b), is shown in figiares 2(a) and 2(b) for the various vortex 
positions discussed. 

, Chordwlse load distribution .- The loading on the tail can be 
calculated from the relation 



2 

YqCL Bx 


2 dui 

YqO. 


The value of 2lA9it /Stilt follows from equations (107); thus 


(108) 


SAqpit ^ 

y Si2-yi2 L " yti2-ai2 


2V0OKI1 + £w 


(109) 


Hence, 




_ 

^^4 L Voajtai(sinh®7i+sin^i) J 


t^CTi 

dx 


( 110 ) 


Equation (UO) when transformed to the ^ plane gives the surface load 
distribution over the tall due to the presence of the two point vortices 
as well as the Inclination of the tail to the free stream. 


Total lift .- The total lift can be obtained by integrating the span 
loading. Thus, if Lt represents the lift oh the tail. 


qa 


2 r‘ 
Voa J 


Acpt dy 


( 111 ) 


Carrying out this integration yields 


Lt 


2T^ 


= - — - (rj^ sihh 7g cos Ug-Si sinh 7j^cos u^^) 

2qa Vqix 


( 112 ) 
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This result ira-s transformed and plotted in figures 3(a) and 3(h) for the 
values discussed. 


Total- drag .- The total drag Is given by the relation 

Dt = ol^ + Ft 


(113) 


where Ft is the suction force on the tail leading edge. As in previous 
steps (see, e.g., equation (92)) the cadculation of Ft depends on the 
evaluation of the function G. In this case G^ is 


a* 


17^-To^ 


2J1 


1 - 


+\i‘ 


2Vo03t 


(U4) 


The final result, in the plane, for the total drag can be written 


qqt,^ 


it(ai2-n2) + i:. 

2it 



Zn 


( a^a-h +( ai2 -hj2 -a ^ ) y ( a^^ -h^s )2-^l^a3ahia 

( at2-hi2 -r^a a^a -h^a -r^a ) y(aQ_2-hj_2-ri2)2+ltej^ah3_a 


(115) 


A plot of the drag is given in figures il-(a) and 4(b) , 

Chord loading .- As before, in the development of equations (97a) 
and (97b^ for example, the chord loading can he calculated either by 
performing the integration /(Ap/(w)t <3y, or by differentiating with 
respect to x the total, lift (Lt/qa) . These two different approaches 
serve to check each other end both lead to the same result, namely. 


/ 


C-) 

\qa-/t 


dy = 4rt0i 


dJx 

dx 


1 - 


Fy slnh 7i cos coj^ 
VQO0xrt( sinh^7i+sin^C0x) 


(116) 


Center of pressure. - Results for the center of pressure Xc.p. 
(■where Xc.p. = -Mt/Lt) are shown in figures ^{a) and 5(b) . 
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CONCLUDING REMARKS 


When an airplane is slender enou^^^ (in the longitudinal sense) or 
is flying close enough to the speed of sound, the mathematical descrip- 
tion of its attendant flow field is greatly simplified - so much so, in 
fact, that the analysis of whole wing-body- tail combinations is feasible. 
This sinplification comes from the fact that the induced velocities in 
each lateral plane are con^jletely independent of the nature of the air- 
plane or flow field behind the reference plane and are affected by 
disturbances ahead only through the presence of free vortices trailing 
dovmstream from the lifting elements. In the case of a tall, these free 
vortices stream back from the wing trailing edge. 

In this report, special wing plan fbrms were studied: special in 
that they produced flat span-loading curves between the wing tips and 
fuselage. For such -vrlngs, the free trailing vortices were concentrated 
entirely in the region directly behind the wing tips. In general, the 
trailing vorticity would be concentrated predominately in this region. 
The behavior of this trailing vortex system is bounded by the behavior 
of two extreme models: a vortex sheet lying everywhere in the plane of 
the tall, and two laterally symmetric point vortices lying in or above 
the plane of the tall. Each of these models was examined. 

One point vortex was placed in the plane of the tail- at a distance 
from the fuselage in the spanwlse direction determined by replacing the 
rtng-span loading curve by a rectangle of the same height. As shown by 
figures 2 through 5^ the resxxlts for this point vortex were not signifi- 
cantly different from those for the vortex sheet. In either case, the 
presence of the trailing wing vortices reduced by about 4o percent the 
effectiveness of the triangular tail surface in producing lift for the 
range of tail spans and body diameters considered. For the same condi- 
tions the tall drag was reduced only l8 percent. 

For the particular locations chosen for the point vortices, it was 
found that both the lift and dreig decreased as the vortices moved closer 
to the tail. On a percentage basis the decrease was roughly the same. 


^^The assunption's underlying slender wing theory are obviously violated 
along lines such as the leading edge, x = mco> and the Mach wave from 
the tralling-edge-fuselage Juncture, x = cq. Along these lines the 
pressure gradient is discontinuous and (Mo^-l) cpxx is not bounded. 
Similar situations appear repeatedly in the linearized analysis of 
aerodynamic flo\^ phenomena and in each case agreement with experi- 
mental results cannot be einticipated. 
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The position of the center of pressure on the triangular tail was 
Insensitive to the presence of the wing vortex system regardless of the 
vortex pattern chosen. In the extreme case, when the point vortices were 
nearest the tail, the location of the tail center of pressure with refer- 
ence to the tail apex as 5 percent forward of the position obtained when 
the wing was absent. 


Ames Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Moffett Field, Calif., Aug. 20, 1951 
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APPENDIX A 

LIST OF IMPORTANT SYMBOLS 


Co 

Cd 

d 

D 

Eo 

E(k,t) 

F 

F(k,i}r) 


horizontal distance from y^=0 plane to vortex 
( See sketch ( z ) . ) 

ai + ihi 

ai - ihi 

local chord 


characteristic chord 
(See sketch (a).) 

section drag coefficient 
section lift coefficient 


qc 


(- 

\ 1C 


section drag force 

drag force ^ 

complete elliptic integral of the second kind 


(m 


^ — dt 


Incomplete elliptic Integral of the second kind 


suction force at leading edge of lifting surface 

incomplete elliptic Integral of the first kind 
dt 



•k2t2 


l-t'= 


dt 



^/(l-k^^)(l-t2) 


Ko 


k 

k» 


vertical distance from z=0 plane to vortex 
(See sketch (z).) 

con 5 )lete elliptic Integral of the first kind 
dt 



Jo y(l-ko2)(i_t2) 
modulus of elliptic integrals 
y l-k2 
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I 

L 


( J 


si‘=-ri‘= 


V(s^^-ro‘^)(s^-t^ 
t(s^-ro^) 

ko in the plane 

2roJ (s^^-ro^)(s^-t^) 

(s^-ro^)(t\ro^) 

k 2 in the plane { 

^ ti Vsi^-ri^ / 

(g^-rpS) A/u^-t^)(s^^ro^ 

( t ^-ro ^) a/ (s^-a^)(s^g^-ro‘^) 

/ ( o i^-r 1^) ( s i^t 1 

k 4 in the li plane I i— ■ ■ -- 

W (Bi2-gl^)(ti2-ri 

section li^t force 
lift force 


% 

2 \ 


m slope of wing leading edge 

(See sketch (a).) 

M pitching mcment, positive vhen tail is forced down 

Mq free-stream Mach nvnnber 


P 

q. 

1 


ro 


s 


So 


S 

t 

to 


local static pressure 


loading coefficient 


n -Pu ^ 

. 1 y 


free-stream dynamic pressure 




radius of body 
(See sketch (a).) 

distance from x axis to wing leading edge 
(See sketch (a).) 

TnaxiTnuTti value of s 
(See sketch (a).) 

wing area 

distance from x axis to -^Ting trailing edge 
(See sketch (a).) 

maximum value of t 
(See sketch (a).) 
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u,v,w 

Au,Av,Aw 

Vr 

Vo 

x,y,z 

^.p.' 

a 


perturbation velocity components in tbe x,j,z directions 
respectively 

jump in velocities across z=0 plane, u^-u^, ^u“^Z 

respectively 

radial component of perturbation velocity in a y,z plane 

free-stream velocity 

Cartesian coordinates 

distance to center of pressure 
(See sketch (w) and fig. 5.) 

angle of attack of airplane 


3 

sink 7^ 
elnh 7g 

r 

n 

e 

Ao(k,V) 

I, 

P 

Po 

a 

9 

A«P 

(A9)t.E. 


a/Mo^I 

Ihil 


ffjLsinUi 



r^^ sin CJ2 


total circulation about wing section 


(A'P)t.e. 


polar angle in y, z plane 


I [KE(k’,ilr) + (E-K)F(k»,i}f)] 

slope of tall leading edge 
(See sketch (x).) 

complex variable (y+iz) 

polar distance in y,z plane ( •/ 

free-stream density 

distance from x axis to leading edge of tail 
(See sketch (x).) 

perturbation velocity potential 

jump in <P across z=0 plane (^-^z) 

value of A*P at trailing edge 
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argument of elliptic integrals 
V(B^^ro4)Cs2-t2Jy2 


f J ^ 

iITq in plane ^ 


2 4-'"2 

Si -ti 


2tA/ro^-y^ 


VCt^+ro^)^ - it-y^t^ 
ijTg in |i plane 


Vri^-: 


71^= 


A/ti^-yi^ 


s(t^+rp^) 

t(s%ro^) 


o2 ^ 2 

s -rp 


2 2 
s +rp 


te in li plane 

s(t^-rp^) 

t(^-rp2) 

^8 in li plane 



ti 

^4 plane [ -f- 


/ 


2 2 
Si -ri 


V^iflSE 

V 81^ -11^ 




20-, 


[ ri^ -ai^-lii^+*\/(ai^+lii^-ri^)^ + 4ri^hi^] 


2ri^ 


Subscripts 

complex plane res\ilting from tbe application of the Joukowski 
transformation to the physical plane 

boundary conditions given by equations (12) 

boundary conditions given by eqiiations (13) 
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a 

Td 

I 

n 

t 

u 

V 

V 


wake 

"body 

lower surface of z=0 plane 
component normal to leading edge 


tail 


upper surface of z=0 plane 

rolled up yortex 

wing 



62 


NACA TN 2554 


APPEKPIX B 

INVERSIOU OF AN IHTE(31AL EQUATION 


The integral equation 


f(Tli) 


2 ^ 

2 n 



Av(Tlg)dTl2 

^ i“'n 2 


(Bl) 


can be inverted by applying operational techniques. Consider the 
operator 

^ Mb 1 

X-Bi 

where h(r] 1) is a function to be chosen later. Operating On both sides 
of equation (b 1 ) yields 





X-BX ^ 



AV(ti 2)Mb 1 ) 
(Bi-Ba)(X-Bi) 


(B3) 


where the order of the integral and differential signs on ’the ri^t-hand 
side of the equation indicates that the x]z integration is to be per- 
formed first . 


The next step is to reverse the order of integration in the double 
integral term in equation (B 3 ). Since an inherent singularity exists in 
the area of integration at the point Bi =" B2 = however, some care 
must be used in order to obtain this reversal. Designating by k(x) the 
difference betrvreen the term taken first with one order of integration and 
then ^rtth another, thus 

E(X)= ArCuajHCu) 

via Ja (Bx-B 2 )(X-Bi) ( B i~B s) ( X-B 1 ) 

(b4) 

R(x) can be evaluated by isolating the singularity and studying the 
difference only in its vicinity.- Hence, 
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R(x) = , 

G->0 


X+e X+e 
X— e X—e 


dTia 


Av(ng)li(Tu) 

(t]i-i12)(x-tu) 


A»+G ^+€ 


- f ^2 f 41117^ 

iL x% <’’1 


Av(iia)h(Tii) 

n2)(x-Tii). 


(B5) 


Using the transformations t]i = X-ex^ and Tja = X-exa, equation (B5) 
can be reduced to the form 


R(x) = Av(x)h(x) 


;/V^T7Srr/V 


dxi 


--i SLi Xi(x2-Xi) 

= 2Av{x)Hx)[' ^ Zn ^ 

Jq H 

= -jt^Av(x)h(x) 


'Ll Xi(x2-xi) -1 


L-Xi 

1+Xl 


(B6) 


Sixbstitute equations (b 6) and (B^4-) into equation (B3) and there 
results 


/’ 


h(t1i)f(3i)dqi _ 1 


X-Tli 


= - ~ -it^Av(x)h(x) + r Av(T)a)dT|2/' 

2^1- 'Ja (ni-n2)(X-ili) 


or 


Av(x) = 


2 

rtli(x) _ 



h(ni)dfn 

(%-Tla)(x-ili) 



X-Ti 


(B7) 


If h(qi) = ^(b-% )(Tii-a) equation (BT) becomes 


Av(x) = — ■ 

Jty(b-x)(x-a) 


1 /'\v(na)a„3 + 

2 Ja Ja X-Tii 


(B8) 

Since, However, equation (b8) contains both Av and its integral, it 
does not represent a unique solution for At. In order to obtain a 
unique solution, some additional condition must be imposed. If this 
condition is that AT(a) = 0, then the proper choice of h(T]i) in equa- 
tion (B7) is 
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■ Mil) = 

V b-li 

which leads immediately to hhe unique result 


Av(l) 



f(T1l)«3.Tll / b-Tli 

J Tli-a 


(B9) 



MCA ™ 2554 


65 


REFERENCES 


1. Jones, R. T.: Properties of Low-Aspect-Ratio Pointed Wings at Speeds 

Below and Above the Speed of Sound. MCA Rep. 835^ 19^6. (For- 
merly MCA TN 1032 ) 

2. Heaslet, Max. A., Lomax, Harvard, and Spreiter, John R.: Linearized 

Compressible-Flow Theory for Sonic Flight Speeds. 

MCA Rep. 956 , 1950 . (Formerly MCA TN l82h) 

3 . Spreiter, John R.: Aerodynamic Properties of Slender Wing-Body 

Combinations at Subsonic, Transonic, and Supersonic Speeds. 

MCA Rep. 962 , 19118 . (Formerly MCA TW's l662 and I 897 ) 

h. Lomax, Harvard, and Heaslet, Max. A. : Linearized Lifting-Siirface 

Theory for Swept-Back Wings with Slender Plan Forms. 

MCA TW 1992 , 1949 . 

5 . Heaslet, Max. A., and Lomax, Harvard: The Application of Green’s 

Theorem to the Solution of Boundary-Value Problems in Linearized 
Supersonic Wing Theory, NACA Rep. 96 I, 1950. (Formerly 
MCA TW 1767 ) 

6. Cohen, Doris: The Theoretical Lift of Flat Swept-Back Wings at 

Supersonic Speeds. MCA TN 1555> 1948. 

7 . Cohen, Doris: Theoretical Loading at Supersonic Speeds of Flat 

Swept-Back Wings with Interacting Trailing and Leading Edges. 

MCA TN 1991 , 194-9. 

8. Broim, Clinton E.: Theoretical Lift and Drag of Thin Triangular 

Wings at Siipersonic Speeds. MCA Rep. 839^ 1946. 

9 . Heuman, Carl; Tables of Complete Elliptic Integrals, Journal of 

Mathematics and Physics, v. 20, 194-1; PP. 127-206, 

10. Spreiter, John R., and Sanks, Alvin H. : The Rolling Up of the Trail- 

ing Vortex Sheet and its Effect on the Downwash Behind Wings. 

Jour, of the Aero. Sci., vol. I 8 , January 1951, PP^ 21-32. 



66 


NACA IN 255^ 



0.100 

.168 .213 

.256 .290 
.332 .361 

.400 .425 
.526 .545 
.644 .661 
.756 .773 


1.083 1.096 


0.316 

.413 

.473 

.526 


1.91 1.92 1.96 2.00 

2.32 2.33 2.37 2.40 
2.83 2.84 2.87 2.91 

3.33 3.34 3.38 3.41 


0.447 

0.500 

0 

b^ 

0 

0 

0 

« 

0 

0 

0.800 

y/mco 

<< 

0 

0 

y/mco 

y/mco 

y/mco 


0.500 

0.600 

0.700 

0 

0 

CO 

ID 

.542 

.593 

.690 

.791 

.883 

.598 

.649 

.746 

.844 

.933 

.648 

.699 

.796 

.890 

.982 

.698 

.748 

.845 

.933 

1.031 

.797 

.845 

.9 hi 

1.028 

1.126 

.896 

.941 

1.033 

1.123 

1.221 

.993 

1.036 

1.128 

1.217 

1.315 

1.090 

1.130 

1.222 

1.312 

1.412 

1.283 

1.323 

1.413 

1.504 

1.606 

1.482 

1.522 

1.610 

1.700 

1.800 

1.676 

1.723 

1.803 

1.893 

1.996 

2.07 

2.12 

2.20 

2.29 

2.39 

2.47 

2.52 

2.60 

2.69 

2.79 

2.98 

3.02 

3.10 

3.19 

3.29 

3.49 

3.52 

3.60 

3.69 

3.78 



























NACA TN 255^ 


67 


TABLE II.- VALUES OF LOCAL CHORD 
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Figure /. - Relations between the wing area, wing span, and 

body diameter. 
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(b) Constant t^/s^ 
Figure /. - Concluded. 
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(a) Vortex sheet. 



(bj Point vortices. 


Figure 3. - Variation of the tiff on the tail. 
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(T^/mc^ 


(a) Vortex sheet. 



(b) Point vortices. 

Figure 4.- Variation of the drag on the tail. 
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Figure 5.- Variation of the center of pressure location 

on the tail. 
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